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ABSTRACT 

We present high signal-to-noise, moderate spectral resolution (R ~ 2000 — 
2500) near-infrared (0.8 — 5.0 /im) spectroscopy of the nearby T Tauri star TW 
Hya. By comparing the spectrum and the equivalent widths of several atomic 
and molecular features with those for stars in the IRTF near-infrared library, we 
revise the spectral type to M2.5V, which is later than usually adopted (K7V). 
This implies a substantially cooler stellar temperature than previously assumed. 
Comparison with various pre-main sequence models suggests that TW Hya is 
only ~ 3 Myr old; much younger than the usually adopted 8 — 10 Myr. Analysis 
of the relative strengths of the H lines seen in the spectrum yields estimates 
for the temperature and density of the emitting region of T e > 7500 K and 
n e ~ 10 12 — 10 13 cm~ 3 . The thickness of the emitting region is 10 2 — 10 4 km 
and the covering fraction is ~ 0.04. Our derived physical parameter values 
agree with the predictions of the magnetospheric accretion scenario. The highest 
signal-to-noise H lines have profiles that indicate multiple emission components. 
We derive an excess spectrum (above that of the M2.5V template) that peaks 
in the H band. Although our derived veiling values, ~ 0.1, agree with previous 
estimates, the excess spectrum does not match that of current models in which 
this flux is generated by an inner optically thin disk. We suggest that the excess 
flux spectrum instead reflects the differences in atmospheric opacity, gravity, and 
age between TW Hya and older, higher gravity field M2.5 dwarfs. 

Subject headings: techniques: spectroscopic - stars: pre-main-sequence - stars: 
individual: TW Hya 
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Introduction 



The optical and near-infrared (NIR) spe ctra of Classical T Tauri stars (CTTSs) are char- 
acteriz ed by strong H emission lines (see e.g., Gullbring et al.|[l998 ; Muzerolle. Hartmann. fc Calvet 
1998bl ). In the currently accepted paradigm for these pre-main sequence objects, these emis- 
sion lines arise from gas being accreted onto the stellar surface from the inner regions of a 
circumstellar disk along the stellar magnetic field lines. In this magnetospheric accretion flow 
model, the inflowing material i s heated and ionized as i t passes through a standing shock 
close to the stellar photosphere (ICalvet fc Gullbringill998l ). Attempts to confirm observation- 
ally the predictions of these models have generally fo cussed on only a few of the strongest H 
lines, usually Ha, Pa/3, and Br7 (e.g. lMuzerolle. Calvet. fc Hartmannll200ll ; lFolha fc Emerson 
200ll ) . Furthermore, much of the theoretical work to date has focussed on modeling individ- 
ual emission line profiles (at high resolution) and the overall spectral energy distributions 
(at low resolution) of CTTSs, with relatively little emphasis on the fluxes and flux ratios 
of a series of emission lines. Only recently have the relative strengths of the higher order 
members of various H series been combined with their stron ger counterparts to estimate the 



densi ty and temperature of the emission regions in CTTSs (IHerczeg et al.l I2004J : iBary et al. 



20081 ) and thereby constrain theoretical models. A complication in carrying out such com- 
parisons of relative line fluxes to model pred ictions arises from the notorious variability of the 
spectra of CTTSs (e.g. lAlencar et al.ll200ll . and references therein), and therefore the need 
to acquire spectra covering all useful lines simultaneously. Consequently, the location and 
physical conditions of the H line emitting regions are not well constrained observationally. 



At a distance of ~ 51 pc (jMamajekll2005l ). TW Hya is one of the nearest T Tauri stars. 
Althou gh it is a member o f a small group of stars known as the TW Hydra association 
fTWA; iKastner et all Il997| ). TW Hya resides surprisingly far from any known molecular 
cloud. Herbigl ( 19781 ) originally identified this object as a possible post-T Tauri star with a 
spectral class of K7 Ve. However, its optical spectrum exhibits a very strong and variable 
Ha emission line a nd emission lines of He I, [ O I], O I, and Ca II, features that suggest 
that it is a CTTS (IRucinski fc Krautterl 119831 ). iRucinskil (119851 ) noted that TW Hya was a 
strong IRAS source, suggesting that it is su r round ed by a dusty accretion disk. This was 
confirmed by lWeintraub. Sandell. fc Duncan! ( 119891 ). who found TW Hya to be a strong sub- 
millimeter source. The ac cretion disk ar o und TW Hya was subseq u ently directly i mage d 
at various wavelengths by Iwilner. et aD (hoooh . krist et al.l fcoooh . iTrilling et all toQlh . 
Weinberger et al.l (120021 ) and iQi et al.l (120041 ) and is seen approximately face-on. Analysis of 
the spectral energy distributio n (SEP) of TW H ya has indicated that the system contains 
a so-called "transitional disk" (jStrom et al.lll989l ) with two separate co mponents: a largely 
dust free, optically thin, hot, inner reg ion, extending from ~ 0.06 - 4 AU ( ICalvet et al.ll2002l ; 



Eisner et al.ll2006t iHughes et al.ll2007| ). and a cold outer region extending to a radius of ~ 200 



-3 - 



AU. TW Hya is actively accreting material from its disk, with an accretion rate of ~ 4 — 20 x 
10~ 10 M yr" 1 jMuzerolle et al.l I2OO0I : Alencar fc Batalhal I2OO2I : iHerczeg et aD l20Q4h . The 
age of th e TWA has been has been estimated from a variety of methods to be between 5 and 
30 Myr fcastner et alill997l: IWebb et al.lll999l: IWeintraub elTall I2OO0I: batalha et aD 12002 



Makarov. Gaume. fc Andrievsky 



2005 



Barrado y Navascuesl2006 ). For a K7V spectral type, 



comparison of the location of TW Hya on t he H-R diagram wi th pre-main sequence isochrones 
gives an age of TW Hya of ~ 10 Myr ( Webb et al. 1999h. although a few substantially 
younger age estimates can be fou nd in the literature (e.g.. iMakarov. Gaume. fc Andrievsky 
2005c iBarrado y Navascuesll2006l ). Sum maries of the prope rties of TW Hya and th e TWA 
have been given in the general reviews by IZuckermanl (120011 ) and IZuckerman fc Songl (12004 ) . 



In this paper we present new high signal-to-noise medium resolution (R ~ 2000 — 2500) 
NIR (0.8 — 5.0/im) spectroscopy of TW Hya and use these data to reassess its spectral 
classification. We derive a new spectral type, which is substantially later than the canonical 
value. The later spectral type implies a much cooler stellar temperature, and hence a much 
younger age than usually adopted. Subtraction of a scaled, matching spectral template 
from our spectrum reveals numerous (primarily H) emission features. Because the spectra 
simultaneously cover both the Pa and Br series, we avoid the ambiguity introduced by 
temporal variability of the line fluxes and we use the relative line fluxes to determine the 
physical conditions of the emitting gas. The derived temperature, density, and thickness of 
the emission region are found to be in good agreement with the predictions of magnetospheric 
accretion models. 



2. Observations and Data Reduction 

We observed TW Hya at the NASA Infrared Telescope Facility (IRTF) on Mauna Kea 
on 2010 Feb 26 and Feb 27 (UT) with SpeX, th e facility near-infrared medium resolution 



cross-dispersed spectrograph (IRayner et al.ll2003l ). Ten individual exposure of TW Hya, each 



lasting 100 s, were obtained using the short-wavelength cross-dispersed (SXD) mode of SpeX 
on 2010 Feb 27 11:26 UT. This mode yields spectra spanning the wavelength range 0.8 -2.4 
jUin divided into 6 spectral orders. Twenty exposures, each lasting 30 s, were obtained in 
the long wavelength cross-dispersed (LXD2.1) mode on 2010 Feb 26 11:20 UT. This mode 
yields spectra covering 2.2— 5.1fim in 6 spectral orders. The observations were acquired in 
"pair mode", in which the object was observed at two separate positions along the 15"- 
long slit. The slit width was set to 0'.'3, which yields a nominal resolving power of 2000 
for the SXD spectra and 2500 for the LXD2.1 spectra. (At the distance of TW Hya of 51 
pc, the SpeX 0'/3 slit spans ±7.7 AU.) The slit was set to the parallactic angle during the 
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observations. The airmass was about 1.74 for both sets of observations. Observations of HD 
92845, an AO V star, used as a "telluric standard" to correct for absorption due to the Earth's 
atmosphere and to flux calibrate the target spectra, were obtained immediately preceding 
the observations of TW Hya. The airmass difference between the observations of the object 
and the standard was 0.07. The seeing was estimated to be ~ 0'.'7 — 0'.'8 and ~ 0'.'6 at 2.2 
jum on 26 Feb and 27 Feb, respectively, and conditions were clear on both nights. A set of 
internal flat fields and arc frames were obtained immediately after the observations of TW 
Hya for flat fielding and wavelength calibration purposes. 



The data were reduced using Spextool (jCushing. Vacca &: Raynerll2004l ). the IDL-based 



package developed for the reduction of SpeX data. The Spextool package performs non- 
linearity corrections, flat fielding, image pair subtraction, aperture definition, optimal extrac- 
tion, and wavelength calibration. The sets of spectra resulting from the individual exposures 
were median combined and then corrected for telluric absorption and flux calibrated using 
the extracted AO V tellur i c stan dard spectra and the technique and software described by 



Vacca. Cushing fc Raynerl (120031 ). The spectra from the individual orders were then spliced 
together by matching the flux levels in the overlapping wavelength regions, and regions of 
poor atmospheric transmission were removed. The final SXD and LXD spectra were merged 
by scaling the SXD spectra by 1.05 in order to match the flux levels of the LXD data in the 
region of wavelength overlap between 2.2 and 2.4 /zm and then computing a weighted aver- 
age. The weights were determined from the signal-to-noise (S/N) values in the two spectra. 
The final 0.8 — 5.0 /j,m spectrum is shown in Fig.[TJ The S/N varies across the spectral range; 
it is of the order of several hundred across the SXD wavelength range, ~ 50 across the 3 — 4 
[im. range and about 10 in the 4 — 5 /im range. 

We computed synthetic NIR magnitudes in various filters from our final spectrum. The 
estimated 2MASS magnitudes are J = 8.35, H = 7.67, and K s = 7.41. Comparison with the 
2MASS point source catalogue indicates that our synthetic magnitudes are systematically 
lower by ~ 12%. Given the variable nature of TW Hya, it is not clear if this difference 
reflects a real change in the source brightness or a systematic error in our absolute flux 
calibration. However our synthetic J — H and H — K s colors agree to within 0.02 mags with 
the 2MASS values, which indicates that our relative flux calibration is very accurate (see 



also lRayner. Cushing fc Vacca! l2009i . where it is demonstrated that the adopted method of 



flux calibration produces spectra with relative fluxes accurate to a few percent ). 
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3. Analysis 



3.1. Spectral Classification 



Although a spectral type as early as K6Ve has been given for TW Hya (ITorres et al. 



2006), the most commonly quoted spectral type for this object is K7Ve following iHerbig 
(119781 ) . However, determinations of the spectral type of T Tauri stars based on optical 



spectra can be biased by the substantia l excess optical continuum produced by the accre- 



tion process (see e.g., iMora et al.l 120011 ) as well as any differential extinction due to dust 



either locally or along the line of sight. In the NIR, however, the contribution from the 
excess continuum is believed to be small for most CTTSs, and this spectral region is much 
less affected by dust than the optical, while still containing numerous diagnostic lines that 
can be used for spectral classification. With the recent publication of th e IRTF Spectral 
Library (IRayner. Cushing fc Vaccall2009l ; ICushing. Rayner. fc Vaccal 120051 ) . it is possible to 
re-evaluate the spectral types for CTTSs by directly comparing their NIR spectra with those 
for a large stellar sample. Therefore, we have re-examined the spectral type for TW Hya 
based on our SpeX spectrum, which was obtained with the same instrumentation and re- 
duced in exactly the same manner as those in the IRTF Spectral Library. We attempted to 
determine the spectral classification of TW Hya using a variety of methods. 

Comparison of the overall shape 0.8 — 5.0 jum spectrum of TW Hya (Fig. [1]) with 
those presented in the IRTF Spectral Library ( IRayner. Cushing fc Vaccall2009[ ) immediately 
suggests that this object has a spectral type between M1V and M4V, and certainly later 
than K7V. In particular, the broad shallow absorption at ~ 0.85 fim, due to TiO bandheads, 
and the prominent continuum 'bumps' in the H and K bands, due to water absorption at the 
edges of the bands at ~ 1.4, 1.9, and 2.7/xm, suggest a spectral type la ter than about M1.5V 



but earlier than about M3.5V (see IRayner. Cushing &: Vaccal 120091 ) . Such a comparison 



should not be affected by extinction along the line of sight to TW Hya because, at a distance 
of only 51 pc and a location well outside of any mole cular cloud, as well as a nearly pole- 
on orientation, TW Hya is believed to have Ay = (IHerczeg et al.l l2004i see §3.2 below). 
Furthermore, the M stars in the IRTF Spectral Library either have zero extinction or their 
spectra have been corrected for extinction. Nevertheless, T Tauri stars are known to exhibit 
"veiling", in which emission from the surrounding accretion disk provides additional flux 
above that of the photosphere in the IR, while the accretion shock contributes additional 
flux in the UV/optical to the total spectrum from the object. As T W Hya is a "transitional 
disk" system, the NIR veiling is believed to be rather low (< 10%; iJohns-Krull fc Valenti 



2001 



see also §3.3 below). Even if the veiling were substantial, the smooth excess flux 
spectrum expected from the disk emission should not produce the broad absorption features 
that characterize the spectra of early M subtypes. 
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We then independently examined various spectral regions, and compared them with 
the corresponding regions of the spectra in the IRTF Spectral Library. We carried out 
this comparison both "by eye" and using a least-squares fitting procedure (to determine 
the best-fit scaling for each spectrum in the IRTF Spectral Library). As can be seen in 
Figs. |2a] and |2bJ both the spectral slope and the depths of various spectral features in the 
0.8 — 1.35 jum region indicate that TW Hya has a spectral type between M1V and M3V. In 
particular, the presence of weak TiO bands at ~ 0.85 and 0.89 jttm, a weak absorption due 
to FeH at 0.87 //m, and the absorption feature due to H 2 at 1.33 jttm suggest a spectral 
type of early M, rather than K. A similar result is found for the K-band region (Fig. l2dj) . 
in which the spectral shape of TW Hya appears to be very similar to that of an M2.5V 
star. For the H-band, however, we find a slight disagreement between the spectral shape, 
which indicates a spectral type of M2V or later, and the strengths of various discrete features 
(particularly the complex of Fe and Si lines between 1.60 and 1.65 /xm), which indicate a 
spectral type of ~ M1V. The x 2 comparison over the wavelength region between 0.95 and 
1.34 /iin (where veiling is expected to be minimal; see below) yields low \ 2 values for the 
M0.5V, M1.5V, M2.5V, and M3V spectra from the Library, with a formal best match with 
the M1.5V spectrum. However, the overall spectrum does not match that of the M0.5V star 
and the shape of the K band spectral region does not match that of the M1.5V star. There 
is very little difference between the M2.5V and M3V spectra. 

We also measured the equivalent widths (EWs) of numerous features in our spectrum 
and compared them with values determined from direct measurements of the spectra of 



K and M stars available in the IRTF Spectral Library ( IRayner. Cushing fe Vaccal 12009 



Cushing. Rayner. fe Vaccal 120051 ) . (Since our spectrum has the same resolution and com- 
parable S/N as those in the IRTF Spectral Library, a comparison with these values should 
provide a better indication of the spectral type than those given by other authors with dif- 
ferent stellar libraries.) The EWs of some of these absorption features (due to Al I, K I, 
Mg I, Na I, and FeH) in the spectrum of TW Hya are presented in Table [TJ A comparison 
between a subset of these values deemed to be the most sensitive diagnostics of spectral type 
and those derived in an identical manner from the K and M stars in IRTF Spectral Library 
is shown in Figure [3] as a function of spectral type. In this figure, the horizontal dashed lines 
denote the range in the EW (±lo") for TW Hya. The results presented in Fig. [31 particularly 
the FeH and K I 1.253/im EWs, clearly suggest a spectral type of M2.5V - M3.5V for TW 
Hya. 

Because of veiling, one must be careful in interpreting the EWs as spectral type diag- 
nostics: lines may appear to be too weak simply because of the additional continuum flux 
from the disk. Therefore, the effect of veiling on a spectral type determination depends on 
the trend of the EW of a feature with spectral type for normal stars. If the EW of a feature 
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decreases with spectral type, veiling will cause a feature to appear weaker and will result in 
a later spectral classification. If the EW increases with spectral type, veiling will result in an 
earlier spectral classification. However, veiling cannot produce lines that are too strong for a 
given spectral type. Furthermore, t he veiling in TW Hya is expected to be negligible in the Y 
and J bands (jEdwards et al.ll2006l ; ISalyk et al.l 120091 ). For these reasons, we relied primarily 
on the diagnostic lines in these bands and on those lines whose strength increases with later 
spectral type. The latter set of lines should yield a limit on the spectral type. Nevertheless, 
even for the H and K bands, the effect of veiling on the spectral type we determine from 
the EWs should be small. The true EW (EWo) can be computed from the observed EW 
(EW obs ) with 

EW = EW obs (l + r) (1) 



where r is the veiling value. iJohns-Krull fc Valentil ( 120011 ) give a K band veiling value for 
TW Hya of = 0.07 ± 0.04. Even a correction of 10% to the observed EW values will 
not yield a spectral type substantially different from what we determine from the observed 
(uncorrected) EWs. Furthermore, such an increase in the EWs will tend to yield later, not 
earlier, spectral types from the most sensitive lines. 

The effects of veiling can be mitigated to a large extent by considering pairs of closely 
spaced diagnostic lines sensitive to spectral class. The ratio of the EWs of two such lines is 
given by 
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We can express the veiling of one line in terms of the other 

r i — r 2 + Ar 

and if the lines are suitably close in wavelength then Ar << r%- Therefore, 
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and the observed EW ratios will reflect the true EWs. If this ratio is sensitive to spectral 
type, it can be used to classify the star. Fig. H] demonstrates the use of this technique for 
TW Hya using the Na I 1.14 /jm and Mg I 1.18/xm lines. This figure demonstrates again 
that the spectral type for TW Hya is approximately M2.5V, independent of any possible 
veiling. The agreement between this determination and that derived from the individual 
lines demonstrates that the veiling around 1 /im must be negligible in TW Hya. 



Based on the results of these three comparison methods, we conclude that TW Hya has a 
spectral type of M2.5V or M3V, with a possible range between Ml. 5 V and M3.5 V. For this 
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spectr al type, the calibrations for M stars given by lLeggett et al.l (119961 ) and lReid fc Hawley 
( 120051 ) yield an effective temperature of 3400±200 K and a mass of O.4±O.1M . To estimate 
the stellar radius, we adopted the 2MASS J b and magnitud e of 8. 22, a distance of 51 pc, 
and the J-band bolometric corrections from the iBaraffe et al.l ( 119981 ) models to compute the 
bolometric luminosity. The radius was then determined from 



R* 



( 



0.5 



With this procedure, we find Lboi = 0.19 ± O.O3L and R* 
luminosity is in g ood agreement with other estimates (0.22 - 
Webb et al.l[l999h . 



(5) 



1.3 ± 0. 2iU The derived 



O.23L ; iKastner et al.lll997 



3.2. Excess Flux - Emission Line Spectrum 



As can be seen in Fig. 1, the NIR spectrum of TW Hya consists of numerous emission 
lines atop a photospheric spectrum. The emission lines presumably arise from the accretion 
process. We generated the spectrum of this emission component by scaling the spectrum of 
an M2.5V star (Gl 381) from the IRTF Spectral Library to that of TW Hya and subtracting 
it. The scale factor was determined to be 0.29 via a least squares calculation between 0.95 and 
1.34 fim. This wavelength region was chosen for two reasons: (1) as stated above, previous 
i nvestigations have found that the veil ing in the Y band (i.e., around 1 /im) is negligible 
( lEdwards et al.ll2006l : ISalyk et al.ll2009l ); (2) matching a library spectrum at any longer or 
shorter wavelengths resulted in substantial regions with negative 'excess' fluxes and veiling 
factors. Since the extinction toward TW Hya is negligible, the latter result is not physically 
justified. The result of the subtraction of the scaled M2.5V template spectrum is shown at 
the bottoms of Figs. I2ati2dl Numerous lines of the hydrogen Paschen (Pa), Brackett (Br), 
Pfund (Pf), and Humphreys (Hu) series, as well as emission lines of O I and Ca II appear in 
the subtracted spectrum. We see n o evidence of the H 2 emission at 2.1218 /im reported by 
Weintraub. Kastner. fc Baryl (120001 ) . The measured line fluxes are presented in Table |2j We 
now discuss the properties of the H, He, and O lines. 



3.2.1. HI 



As shown by iBary et al.l ( 120081 ) for a sample of 15 classical T Tauri stars and by 
Najita et al.l ( 120101 ) for the case of TW Hya itself, the relative H line strengths in the spec- 
tra of these objects can be well reproduced with the Case B recombination approximation. 
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Although the relative line strengths within any given series can be used to estimate the 
temperature and density of the emitting re gion, a more sens itive diagnostic is provided by 
the inter-series ratios of line strengths (e.g., iBary et al.ll2008l ). Therefore, we computed the 
flux of each line relative to that of Pa 0, the strongest li ne in our spectrum . We then 



used the theoretica l H lin e fluxes for Case B as provided by IStorey fc Hummer! ( 119951 ) and 
Hummer fc Storey! (119871 ) to compute the corresponding line ratios for electron tempera- 
tures between 500 and 30000 K and electron densities between 10 2 and 10 14 cm~ 3 . For each 
temperature and density pair we computed a x 2 value between the theoretical and observed 
line ratios. Contours of the \ 2 values are shown in Fig. |5j The minimum y 2 was found for 
T e = 20000K and n e = 10 13 cm -3 and the corresponding theoretical line ratios for these 
values are plotted along with the observed line flux ratios as a function of upper level quan- 
tum number in Fig O However, as seen in Fig. the electron temperature is not tightly 
constrained by the line ratios and we find that any temperature T e > 7500 K can provide 
an adequate representation of the data points. For temperatures 7500 < T e < 12500 K, the 
best fit electron density is n e = 10 12 cm -3 , while for temperatures T e > 12500 K, the best fit 
yields n e = 10 13 cm -3 . In our analysis below we will use both n e = 10 12 cm -3 and n e = 10 13 
cm -3 , with the corresponding temperature values. For comparison, we also plot i n Fig ffl 



the theoretical values for T e = 3000 K and n e = 10 10 cm 3 deri yed by IBary et al.l ( 120081 ). 



This comparison clearly demonstrates that the values found by IBary et al. J2008I ) cannot 
reproduce the observed line ratios for TW Hya. W e also analyzed the H line strengths seen 
in the Spitzer IRS spectra of TW Hya presented by lNajita et al.l (120101 ) in the same manner. 
Because of the large error bars and relatively few data points, the H flux ratios relative to 
Hua in the SH2 spectrum provide only weak constraints on the electron temperature and 
density. Although the formal best fit for the SH2 line fluxes is found for T e = 30000 K and 
n e = 10 9 cm" 3 , there is very little difference between the theoretical line flux ratios for this 
se t of paramete r s and the values we found from our NIR line analysis. For the SHI spectrum 
of iNajita et al.l (120101 ). t he best fit was fou nd for T e = 5000 K and n e 
the values advocated by IBary et al.l ( 120081 ). 



10 10 cm 3 close to 



As stated above, it is the combination of multiple H line series that provide the strongest 
constraints on the temperature and density. With just the Pa H series alone, for example, 
it would be nearly impossible to distinguish lines produced in a region with n e = 10 4 cm -3 
from those produced in a region with n e = 10 13 cm" 3 . However, the fluxes of the H lines 
in the spectra of T Tauri stars are known to vary, and our LXD data were not acquired 
simultaneously with the SXD data (the time difference was almost exactly 24 hours). There- 
fore, we also carried out our analysis using only the Pa and Br lines in our SXD 0.8-2.5 jum 
spectrum. We obtained results similar to those given above: a temperature of 3 x 10 4 K and 
a density ~ 10 14 cm -3 provided the best fit in this case, but T e = 20000 K, and n e = 10 13 
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cm -3 were nearly equally acceptable. This, along with the fact that the observed flux levels 
of the SXD and LXD spectra differed by only 5% between the two nights, suggests that TW 
Hya did not vary substantially between the two observations, and confirms the results found 
by including the Pf lines in the analy s is. W e note that the H line strengths in the spectrum 
of TW Hya obtained by lCovey et al.l (120101 ) in 2008 are only about 15-20% larger than the 
values we measure. 

The measured H lines contain several sets for which the N upper is the same (e.g., Pa (3 
and Br a, or Pf (3 and Br 7 and Pa 5). The flux ratios of these lines are fairly insensitive to 
the density and temperature of the emitting gas, particularly at densities above 10 9 cm -3 , 
and hence can be used to estimate the extinction. The agreement with the Case B ratios 
shown in Fig. [6] indicates that the assumption of zero reddening, and hence low dust content 
in the emitting region, is correct. 

We us ed the flux in the Pa and Br 7 l i nes to estimate the accretion rate using the cali- 
brations of lMuzerolle. Hartmann. fc Calvetl (jl998bl ). These lines yield accretion luminosities 
of L acc = 5.9 x 10~ 2 L Q (from Pa 0) and 2.9 x 1O _2 L (from Br 7), and mass accretion rates 
of 6.1 and 2.9 x lO~ 9 M yr _1 , respectively, from the relation 



M 



~gm7 



which can be expressed as 

logM [M yr" 1 ] 



log 



£0 



+ log 



R Q 



- log( 



M 



- 7.50 



(6) 



(7) 



These values are considerably higher than that es timated by iMuzerolle et al. 



TW Hya, but are well within the range found by iNajita. Strom, fc Muzerolld (120071 ) for 



(2000) for 



transitional disk systems in Taurus. Furthermore, the ac cretion luminosities are well within 
the ranges of values determined by iBatalha et al.l (120021 ) from their analysis of optical data 
for TW Hya. Under the assumption that the Pa ft line is optically thin, the line luminosity 
also provides an estimate of the ionization rate from 



Qo 



e Pa/3 



(8) 



or 



Qo = 3.8 x 10 



13 



J Pa/3 



10- 13 / vio- 26 



CPa^ 



ph s 



(9) 

We find the ionization rate to be Qo ~ 4.5 — 6.6 x 10 42 ph s _1 for n e = 10 12 cm -3 and 



Qo ~ 6.1 - 7.0 x 10 42 ph s" 1 for n 



10 13 cm 3 . 



The values derived for a density of 10 



12 



- 11 - 



cm 3 are in good agreement with that derived by lNajita et al.l (120101 ). while the values found 



for n e = 10 cm are about a factor of two larger. 
The volume of the ionized region is given by 

Qo Lp a p 



V 



nla B 



epa/^ 



or 



V = 3.8 x 10 



35 



io- 26 



n,. 



10 



12 



cm 



(10) 



11) 



where we have assumed n p ~ n e . The emitting volume was found to be surprisingly small, 



V ~ 0.6 - 4.0 x 10 31 cm 3 for n P 



10 12 cm" 3 and 1.4 



3.3 x 10 29 cm 3 for n P 



10 13 cm 3 . 
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We derive a mass of the emitting gas, mnTipV, of between 1.1 — 6.2 x 10 19 g (for n e = 10 
cm -3 ) and 2.3 — 5.6 x 10 18 g (for n e = 10 13 cm" 3 ). These va lues are substantially lower 
(factors of 10 4 — 10 5 ) than that derived by iNajita et al.l (120101 ) . who assumed a density of 



n, 



< 10 s cm 3 . Assuming the emitting volume is located at or near the surface of the star, 



we can estimate the thickness of the emitting region from 

AR ] 



(12) 



where R* is the radius of the star and /# is the fraction of the surface area covered by the 
accretion spot. Combining this with Eqn. [TTJ we have 



AR = 6.3 x 10 £ 



/ R* 
\R^ 



-2 



CPa/3 

To 



-26 



-1 



n. 



10 



12 



0.01 



km 



(13) 



If we require that the optical depth in the Pa (3 line rp a/ g < 1, as suggested by the good fit 
of the theoretical Case B values to the obs erved line ratios, we can p lace a lower limit on 
Using the line opacity factors, Q, given by lStorey fc Hummer! (119951 ). we have 



/* > 



n 



Pa^ -kPa/3 



e Pa/3 AttRI 

and we find /* > 0.01 for n e = 10 12 cm -3 and /* > 0.02 fo r n e 



10 13 cm 3 . 



(14) 



These limits are 



in reasonable agreement with the values of /* derived by iBatalha et al.l (120021 ) . Using their 
relation between /* and L acc 



/, = O.49(L acc /L ) + 0.01 



(15) 



we derive /* = 0.04 from the Pa (3 luminosity. Using this value in Eqns. [12] and [131 we derive 
the thicknesses corresponding to the emission volumes of 1.6 — 9.1 x 10 4 km for n e = 10 12 



cm 3 and 3.4 — 8.2 x 10 2 km for n P 



10 13 cm 3 . 



The ranges of thicknesses in turn imply 
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column densities of N H ~ 1.6 — 9.1 x 10 21 and 3.4 — 8.2 x 10 20 cm -2 , respectively. Since TW 
Hya is seen nearly pole-on, the relatively small volumes and thicknesses we derive suggests 
that the H lines arise primarily from a fairly thin region, consisting of one or more hot spots 
located presumably at the base of the accretion column at the surface of the star near the 
magnetic pole. 

In Fig. [7] we show the profiles of the Pa and Br 7 emission lines. Overplotted is 
a Gaussian with a FWHM of 150 km s— 1. (Recall that our spectra have R = 2000 in 
this wavelength range, which corresponds to a velocity resolution of ~ 150 km s -1 .) The 
high S/N afforded by our spectra allows us to discern multiple components in these profiles, 
which are replicated in many of the stronger H lines: a narrow (unresolved) core at the rest 
wavelength, a weaker component that extends over a fairly large velocity range (from roughly 
—350 to +250 km s _1 ), and an inverse P Cygni absorption feature on the red side at ~ 300 
km s -1 . We note that extended wings are not seen in the (substantially weaker) O I and 
Ca II line profiles, whereas they are present on even the weak, higher order Paschen lines, 
a result that suggests that they are not due to our point s pread function. In addition, the 
Pa 5 line was used to generate the telluric correction (see IVacca. Cushing fc Raynerl 120031 ) 
and so should be the least affected by any residuals from the telluric correction process, yet 
it clearly exhibits broad wings. Also, the broad emission wings are not seen on many of 
the H lines in the reduced spectrum of TW Hya, but appear only after subtraction of the 
M2.5V template. Furthermore, b road emission over a range of velocities up to —400 km s" 1 



is not uncommon in CTTSs (e.g Edwards et al.l 120061 ) . The H line profiles suggest that the 



strongest emission arises at the base of the accretion column (zero velocity), with additional 
emission from an outflowing wind with velocities up to about 300 km s , and absorption 
from gas that is falling onto the star with a velocity of ~ 300 km s -1 . 



3.2.2. He I 



Exhibiting two absorption troughs on either side of an emission peak at the line center, 
the He I line profile (Fig. [8]) is more complicated than those of the H lines. The blue shifted 
absorption component has a minimum at ~ —100 km s -1 and extends to about —300 km 
s _1 , while the weaker redshifted component has a minimum at ~ 300 km s _1 and extends to 
about +500 km s -1 . The velocity of the red absorption feature is the same as that seen in the 
H line profiles, which suggests that the H and He absorption features are formed in the same 
region of the accretio n flow or infall onto th e star. Our He 1 1.083 fim line profile is unlike any 
of those presented by Edwards et al.l (120031 ). who observed the He 1 1.083 /xm li ne in a sample 
of six YSOs. Nor is it similar to the He I 1.083 /im line profile presented by iDupree et al. 
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( 120051 ) from their observations of TW Hya in 2002. In the Edwards et al.l ( 120061 ) sample 
of 39 CTTSs (which includes TW Hya), the only profile the current observations of He I 
1.083 /im resembles is that of CY Tau, although the extent in velocity of the blueshifted 
absorption we observe (with our limited velocity resolution of ~ 150 km s " 1 ) agrees with 



that o f the blueshifted emission seen in the TW Hya spectrum presented by Edwards et al. 
( 120061 ). These comparisons indicate that the He I line profile in TW Hya is remarkably 
variable. This conclusion is strengthened b y the profile of He I 1.083 /im in the recently 
published spectrum of TW Hya obtained by I Covey et al.l (120101 ) in 2008 using SpeX at the 
IRTF. The He I 1.083 /im profile in their spectrum exhibits a considerably stronger emission 
peak, a somewhat shallower P Cygni absorption trough, and almost no evidence of the red 
absorption feat ure. This change in the He I line profile suggests a denser outflow during the 
observations of ICovey et al.l (120101 ). which would be consistent with their somewhat stronger 
H lines. 



3.2.3. OI 



As discussed by iGrandil (119751 ) and summarized by iRudy. Rossano. fc Puetterl (119891 ) , 
the permitted O I lines seen in the spectrum of TW Hya can arise from recombination, 
fluorescence excited by the stellar continuum, or resonance fluorescence excited by Ly 0. 
Because continuum fluorescence should generate an 0.8446 /im line that is considerably 
stronger than the 1.1287 /im line, while the observations indicate that these two lines have 
approximately equal strength, this excitation mechanism can be ruled out. The relative line 
strengths produced by recombination can be derived from the volume emission coefficients, 



J8446 — n e n Oir 



8446 eff 



4tT 



a: 



8446 



Jii287 = n e n on - 



1287 eff 



47T 



a 



11287 



and 



(16) 
(17) 

4?r -Pa/3 • ( 18 ) 

Because O I and H I have very similar ionization potentials, we might expect O and H to 
have similar spatial distributions and ionization fractions in the emission regions. In this 
case, the strength of the 0.8446 /im line relative to Pa /3 can be expressed as 



JPa/3 = n e n H n- 



hvp, 



a/9 eff 



-a 



J8446 



n O \ Ap a ff«8446 



eff 



(19) 



JPa/3 V ^H^A 8 446a P 2 

This is provides an upper limit to the flux ratio because it assumes all the O is in the form 
of O II; some O may be in O III and higher ionization stages. 
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Extrapolating the effective recombination coefficients of lEscalante fc Victor! (11992! ) for 
O I to electron t e mper atures of ~ 10 4 K, and using the H line emissivi ties given by 



Storey fe Hummer! (119951 ) and a solar oxygen abundance (lAsplund et al.ll2009l ). we estimate 



line strength ratios resulting from recombination of j8446/Jii287 ~ 2 and jsue/jp&p ^ 2.6 x 
1CT 3 . The observed line ratios are found to be F 8 44 6 /Fn 2 87 pa 1.4 and F 8 4 46 /Fp a/ 3 pa 7.7 xl0~ 2 . 
Although it could be argued that the former ratio could be affected by reddening, which could 
be responsible for reducing the observed ratio below the expected value, the latter ratio is 
already ~ 30 times stronger than predicted and would only increase with any corrections for 
reddening. Furthermore, to bring the former line ratio into agreement with the predictions 
would require a reddening correction corresponding to Ay pa 1.9, which seems far too large 
for TW Hya. Therefore, under the assumption that the H and O I emission lines are gen- 
erated in the same region, we can rule out recombination as the source of the O I emission. 
This leaves fluorescence by Ly /3 as the excitation mechanism for the O I lines. This mecha- 
nism should produce equal numbers of photons in the 0.8446 and 1.1287 /im lines, which is 
exactly what we observe. The fact that we do observe equal numbers of photons in the lines 
again confirms that there is no differential reddening along the line of sight to the emitting 
region and yery little dust within the emission region itself, as found previously by many 
others (e.g. Muzerolle et al.ll2000l ). 



3.3. Excess Flux - Continuum Spectrum 

In Fig. [9] we present the 0.8 — 5 /im spectrum of TW Hya with the scaled IRTF Library 
spectra of a K7V star (HD 237903) and an M2.5V star (Gl 381) overplotted. The middle 
panels show the result of subtracting the scaled library spectra from the observed spectrum, 
while the bottom panels present the ratios of this excess spectrum to the scaled library 
spectral templates. Both library spectra were scaled to match the observed spectrum of 
TW Hya between 0.95 and 1.34 /im. Template spectral types of M2V - M3.5V yield very 
similar excess emission spectra to that shown in the middle right panel of Fig. As seen 
in the middle left of Fig. [HI earlier template spectral types yield excess emission spectra 
with several broad "features" , particularly between the J, H, and K photometric bands, and 
somewhat stronger excess emission longwards of ~ 3 /im, but wholly unphysical negative 
values shortward of ~ 1.1/im. 

The excess flux spectrum we derive from subtracting the M2.5V star exhibits a rise to 
shorter wavelengths (A < 1 /im) that is most likely due to the (possibly optically thick) 
continuum emission produced by the accretion process ("optical veiling"). The excess at 
longer wavelengths could be due to the thermal emission from the accretion disk, either 
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the optically 



( lEisner et al 



thick disk at R > 4 AU (ICalvet et al.l 120021 ) or an optically thin inner disk 



2006) . However, th e excess spectrum cannot be fit with any simple disk model; 



the model of lEisner et al.l ( 120061 ) is shown as a dotted line in this panel and clearly does not 
match the shape of the derived spectrum, greatly underestimating the flux in the H and 
K bands and overestimating the flux at longer wavelengths. The peak in the excess flux 
spectrum in the H band requires a dust temperature of ~ 1600 K, but optically thin dust 
disk models, with power law distributions for the temperature and surface density, generally 
cannot reproduce the sharp drop off in flux on either side of this peak that is seen in our 
data. 

An 'excess flux' spectrum that peaks in the H band is very reminiscent of the spectral dif- 
ference observed between late type giants and dwarfs. As can be seen in t he series of spectra 
presented in the IRTF Spectral Library (IRayner. Cushing &: Vaccal 120091 ) . for a given spec- 
tral type a late type giant star exhibits a prominent H band excess over that of a dwarf star 



due to the decrease in the H opacity at the lower density and gravity (jWing fc J0rgensen 



20031 ; iMihalaa Il978l ) . To demonstrate this, in Fig. [10] we plot the 'excess' spectrum derived 
by scaling the spectrum of an M2.5V star to that of an M2.5III star, both taken from the 
IRTF Spectral Library, over the wavelength range 0.9 — 1.34/im, and subtracting the for- 
mer from the latter. The resulting difference spectrum was then scaled to match the level 
of the H band excess in the spectrum of TW Hya and overplotted on the TW Hya excess 
spectrum. As can be seen from this figure, the shapes of the excesses are nearly identical, 
with the giant — dwarf excess extending from the H band through the K band and into the 
M band. Clearly, of course, TW Hya is not an M giant. Nevertheless, the mass and radius 
values determined above yield a gravity of log g = 3.8, far below the valu e of ~ 4.8 estimated 



for fie ld M2.5 dwarfs from the parameters given by lLeggett et al.l (Il996[ ) and lReid fc Hawley 
( 120051 ) . Therefore, we believe the TW Hya H and K band excess is in fact not indicative 
of disk emission, but rather is a signature of TW Hya's youth (relative to field M dwarfs), 
as it contracts and evolves down to the main sequence. This conjecture might be tested 
with higher resolution observations of atomic absorption li nes, whose widths would reflect 
the va lue of logg, in much the same manner as described by lCruz. Kirkpatrick. fc Burgasser 
( 120091 ) for L dwarfs. If the excess could be calibrated as a function of log g, it could provide 
a separate estimate of the age of the system. 

Under the assumption that the excess spectrum is due to emission arising from the 
accretion process, the ratio of the excess spectrum to the library template yields the veiling 
r. As seen in the bottom right panel of Fig. [HI the veiling derived from the M2.5V spectral 
template is fairly small: ~ 6% between 0.8 and 0.9 fim, ~ 8% between 1.4 and 1.8 /im, ~ 10% 
between 2.0 and 2.4 /im, and ~ 12% bet ween 3.0 and 4.0 /xm. Our deri ved veiling values are 
in good agreement with the estimate of iJohns-Krull fc Valentil ( 1200 ll ) for the K band and 
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those of lAlencar &: Batalhal (120021 ) for 0.8 — 0.9 /mi. However, if the suggestion given above 
for the shape of the excess spectrum is correct, only a fraction of this excess emission arises 
from the accretion process and therefore our veiling estimates for the H, K and L bands are 
strict upper limits. Although it is uncertain how much of the excess longwards of ~ 1 /jm is 
due to the effects of the lower gravity of TW Hya relative to the comparison field M dwarf 
used for the subtraction, it is possible that this difference in gravity is responsible for nearly 
all of the observed excess in the H and K bands, and a significant fraction of the excess in 
the L band. 



4. Discussion 



4.1. Spectral Type 



The wavelength region around 1 fim is ideal for spectral classification of CTTSs: (1) the 
excess emission arising from the accretion process, which is very strong in the UV and optical 
wavelengths, is relatively small; (2) the emission from the accretion disk itself, which can 
dominate the spectrum at K band and longer wavelengths, is minimal; (3) the effects of any 
extinction are mitigated substantially compared to the optical; and (4) there are numerous 
strong photospheric absorption features present that can be used for classification. The last 
point is particularly important because, while the contributions from the accretion shock 
and circumstellar disk (UV and IR excesses, respectively) are at a relative minimum in the 
Y and J bands, they are not neces sarily negligible and evidence for significan t J band veiling 
has been found in some sources jFolha fc Emerson! Il999l : ICieza et al.ll2005f ). However, the 
ratio of absorption EWs, as discussed in §3.1, provides a veiling independent method of 
determining the spectral type of a source. Although the relatively 'cle an' nature of this 
spectral region has been pointed out before (IKenyon fc Hartmannl Il990l) and exploited to 



study the emission associated with the accretion process in CTTSs (jEdwards et al.l 12003 



20061 ). this is the first time of which we are aware that the absorption features in this region 
have been used for spectral classification of a CTTS. The absorption line EWs we measure, 
in comparison with those derived from the IRTF library spectra, provide strong evidence 
that the underlying star in TW Hya has a spectral type of ~ M2.5V. This result suggests 
that the values for the stellar mass an d radius typically adopted in the literatur e (T pff ~ 4000 
K, M ~ O.6M , and R ~ 1R Q ; e.g., IWebb et al.lll999l : iMuzerolle et~al] bood feisner et al. 



20061 ) are incorrect, as they are based on an incorrect spectral type (K7V). We suggest that 



the appropriate values for TW Hya are T eff ~ 3400 K, M ~ 0.4M Q , and R ~ 1.3-R^, based 



on an M2. 5V spectral type , the t emperature and mass calibrations given by iLeggett et al. 
( 119961 ) and iReid fc Hawleyl (120051 ) . and the estimated bolometric luminosity. 
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The spectral type we determine, and the stellar parameters we derive, are clearly in- 



consi stent with those found by others (e.g. iWebb et al 
20051 ). Examination of the spectra shown by 



1999 



Yang. Johns-Krull. fc Valenti 



Webb et all (Il999h indicates that the K7V 



spectral classificatio n is appropriate, as TW H ya c learly has an earlier type tha n thos e of 
the M stars shown, klencar & Batalhal J2OO2I ) and lYang. Johns-Krull. k Valentil f|2005h fit 
atmosphere models to high resolution optical spectra and derived parameters again con- 
sistent with a K7V star. On the o ther hand, the mass estimated for TW Hya from the 

Yang. Johns-Krull. fc Valentil (120051 ) stellar parameters is far too large for eve n a K7V spec- 

tral ty pe. Furthermore, the moderate resolution optical spectrum presented by lRucinski &: Krautter 
(119831 ) (their Figure 9), with its strong rise in the continuum to the red and its sharp and 
prominent TiO bandheads near 7100 A, is more suggestive of an early M dwarf than a late K 
d warf (see below) and does not seem to be completely consistent with the spectrum presented 
by lWebb et al.l ( 119991 ). Can these two discrepant spectral type classifications be reconciled? 
The optical spectrum of TW Hya suffers from considerable and variable veiling (see e.g., 
Alencar fc Batalhal 120021 ). so it is possible that veiling has contributed to mis-classification, 
especially if the spectral comparison is restricted to very small wavelength regions lacking 
strong lines diagnostic of spectral type. To investigate thi s possibility, we selected an M2V 
star from the Indo-US Spectral Library ( jValdes et al.ll2004l ) and added veiling in an attempt 
to reproduce a K7V star s p ectru m between 6400 and 7600 A (the spectral region displayed 
in Fig. lb of IWebb et al.l (I1999I )). We found that it was n ot possible unless the veil ing 
was substantially larger (r ~ 3 ) tha n that estimated by either lAlencar fc Batalhal (120021 ) or 



Yang. Johns-Krull. fc Valentil (120051 ) in this wavelength range. Therefore, veiling cannot be 
completely responsible for the classification differences. We also considered the possibility 
that TW Hya is an unresolved binary. However, given the relative luminosities of K and M 
stars, it is not physically possible for a coeval K7V+M2.5V binary system to generate a spec- 
trum in which the K7V component dominates the optical spectrum but is nearly invisible in 
the NIR spectrum. 

Because both the continuum and line emission are know to vary substantially in TW 
Hya, it is tempting to consider the possibility that the spectrum (and spectral type) of the 
object varies as a function of both accretion rate and wavelength. A physically plausible 
explanation for this variation is that the spectrum of TW Hya consists of two separate com- 
ponents, one due to the star (M2.5V) and one arising from an accretion hot spot, either 
on the surface of the star or located in the inner disk, or from the inner disk itself. In 
this scenario, the hot spot or inner disk would have a temperature of ~ 4000 K and would 
generate a spectrum similar to that of a K7V star. This component would also vary in 
strength with the accretion rate, dominating the (optical) spectrum and completely over- 
whelming the underlying M2.5V stellar spectrum when the accretion rate is high. Support 
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for this suggestion can be found by comparing the observed optical spectra with the reported 
Ha equivalent wid ths, a proxy for the accretion rate. The optical spectrum presented by 
Webb et all Jl999h is clearly th at of a K7V star and yields an EW(H a) = -220 A. Similarly, 



the optical spectrum shown by lYang. Johns-Krull. fc Valentil ( 120051 ) is consistent with that 
of a K7V star and has EW(Hoj ) ~ —146 A (C. M. Johns-K rull, private communication). The 
optical spectrum presented by lRucinski fc Krautterl ( 119831 ). however, corresponds to spectral 
type intermediate between a K7V and an early M star, and has EW(Ha) between ~ —71 
and —101 A. In fact, when we attempted to reproduce the features seen in this spectrum we 
found that a combination of a K7V and an M2.5V, with roughly equal contributions at 0.8 
/xm, provides a much better representation than any single K or early M spectrum alone. 
The shape and depth of the TiO bands at 7100 A and the simultaneous presence of these 
bands and a strong Ca I line at 7148 A indicates that the observed spectrum is a composite 
of both a K7V and an early M type spectrum. In our NIR spectrum of TW Hya, the region 
shortward of 1/jm can also be well fit with an M2.5V spectrum and a small contribution 
from a K7V spectrum, comprising ~ 30% of the total flux at 0.8 /im. This is shown in Fig. 
[TlT where it can be seen that at wavelengths shorter than ~ 0.6/im we would expect that, 
during our observations of TW Hya, the K7V spectrum contributed roughly half of the total 
flux, but comprised only a very small fraction of the flux longward of 1/im. Unfortunately a 
more satisfactory test of this suggestion cannot be carried out with the spectra and accretion 
rates or EW values found in the literature. (Surprisingly few optical and NIR spectra of TW 
Hya are actually shown in the voluminous literature on this object.) Frequent simultaneous 
monitoring of TW Hya in both the optical and NIR wavelength regimes would be able to 
confirm or rule out our suggestion. 



4.2. Age 



Our results for the spectral type and effective temperature require a re-evaluation of the 
age of TW Hya. The age adopted for TW Hya under the assumption of a K7V spectral type 
with T e ff = 4000 K is often ~ 10 Myr, d erived by comparin g the location of TW Hya in the 
HR diagram with model isochrones (e.g. I Webb et al.lll999l ). In Fig. [T2J we plot the location 
of TW Hya on the H-R diagram, using our estimated T P ff , along with the model tracks and 
i sochr ones of Baraff e et all (119981). iD'Antona fc Mazzitellil (119981) , ISiess. Dufour. fc Forestini 
(hoool ). and lTognelli. Prada Moroni, fc DegFInnocentil (hold ). Despite the considerable vari- 
ation in the model predictions, all of the models yield masses smaller than our adopted value 
of 0.4M Q and ages substantially younger than the canonical value of 10 Myr. Large under- 
estimates of both mass a nd age by the various models have been pointed out previously by 
([Hillenbrand et all 120081 ) . However, even accounting for an underestimate of the age by a 
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factor perhaps as large as two ( Hillenbrand et al.l 120081) would be insufficient to make the 
star 10 Myr old. The models of iBaraffe et al.l (119981 ) provide a mass of 0.35 < M < O.4M , 
closest to our adopted value of O.4M . For these models, the age of TW Hya estimated 
from the isochrones is ~ 3 Myr. We note that the new location of TW Hya on the H-R 
diagram actually brings this so urce into bet t er agr eement with other members of the TWA, 
as can be seen from Fig. 3 of IWebb et al.l ( 119991 ). which indicates that, aside from TWA 
6 and TWA 9 A and B, the locations of the members of the TWA on the H-R diagram 
are consistent with an isochrone age for the association of ~ 3 — 5 Myr. This age is also 
in reasonable agreement with th e expa nsion age of the TWA of 4.7 ± 0.6 Myr derived by 



Makarov. Gaume. fc Andrievskyl ( 120051 ). and therefore alleviates to a considerable degree 



the previous discrepancy between the isochrone age of TWA and dynamical age estimates. 
Furthermore, compared to the usually adopted age of 10 Myr, an age of 3-5 Myr for TW 
Hya is in much better agree ment with estimate s of the half life of accretion disks in low mass 
stellar systems of < 3 Myr (jEvans et al.ll2009l ). By 10 Myr the disks in such system s have 
usually transformed into debris disks with masses of ~ 0.01 — 0.1M® (jWyattl 120081 ). The 
large mass of the disk in TW Hya (0_.0_6Mq ICalvet et al.ll2002l ) therefore suggests an age of 
considerably less than 10 Myr (jWyattll2008l ). in agreement with our much younger estimate. 
We also note that TW Hya exhibits many of the same emis sion features, and ha s exactly the 
same Li 6708A absorption EW, as the M2.5 star GN Tau (TWhite fc Basra l2003| ). The latter 
star is a low-mass T Tauri star with an estimated age of ~ 3 Myr, again providing support 
for our age estimate for TW Hya. 



4.3. Accretion Parameters 



The determination of an accurate spectral type permits the subtraction of the correct 
spectral template to reveal the emission from the accreting gas and the circumstellar disk. 
Under the assumption that the H lines are optically thin and follow the Case B predictions, 
our analysis of the H line flux ratios yields a temperature of ~ 20000 K, a density of 10 13 
cm -3 and a thickness of the emission column of AR ~ 10 3 km. Comparison of our results 



with those presented by ICalvet fc Gullbringl (119981 ) . who modeled the emission from the 
accreting column of gas in the magnetospheric infall scenario, reveals surprisingly good 
agreement for some of the physical parameters. The models of ICalvet fc Gullbringl (119981 ) 
include emission from the shock and post-shock region, the pre-shock region, and the heated 
stellar atmosphere. The emission properties of the models are parametrized in terms of the 
energy flux carried by the accretion column into the shock, J 7 , given by 



J 7 = 9.8 x 10 



10 



M 



M 



lO-SMeyr- 1 / VO.5M / \2R Cl 



R 



0.01 



erg cm 2 s . 



(20) 
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The H line emission in these models arises from the pre-shock region, which is predicted to 
have a density given by 



n H = 5.8 x 10 



12 



M 



lO-SMoyr- 1 / V0.5M G 



M 



-0.5 



R N-i-V f, 



2R Q 



0.01 



cm 



(21) 



For the stellar parameters for TW Hya given above, these equations give J 7 ~ 4.4 x 10 10 
erg cm -2 s _1 and nn ~ 1.9 x 10 12 cm -3 . The latter value is fairly close to our empirical 
results (under the assumption that all of this material becomes fully ionized in the pre-shock 
region). The models also predict a temperature for this region of T ~ 20000 K, and thickness 
of the ful ly ionized zone of A R ~ 10 3 — 10 4 km, again in good agreement with our findings. 
Recently Eisner et al.l ( 120101 ) has used interferometry to determine the stellocentric radii of 
the Br 7 emitting gas in T Tauri and Herbig Ae/Be stars. They find that the emission 
in general arises from gas within 0.01 AU (1.5 x 10 6 km), an upper limit tha t our results 
certainly satisfy. The infall velocity in the models of lCalvet fc Gullbringi (119981 ) is predicted 
to be 



v 



275 



M 



0.5 



R 



-0.5 



km s 



-1 



(22) 



.O.5M / \2R & , 

or ~ 300 km s _1 for the mass and radius we derive for TW Hya, in excellent agreement with 
the velocities corresponding to the observed inverse P Cygni absorption features seen in our 
H and He I 1.083/im line profiles. 

Using our derived density and placing a limit on the optical depth of the Pa line, we 
derive a lower limit to the fractional area covered by the accretion column of /* > 0.02. Our 
limit is well within the range of previous estimates of the covering; fraction, for both TW Hya 



itself as well as other CTTSs (IHerbst k Koretlll988l : iBatalha et al.ll2002l ; ICalvet fc Gullbring 



9981) . Although our limit is a factor of 7 larger than that estimated by iMuzerolle et al.l 



(120001 ). we also derive substa ntially larger (fac t ors of 4-5) values of the accretion luminosity 



and rate than these authors. (IMuzerolle et al.l (120001 ) also derive a substantially larger value 
of the accretion energy flux J 7 than we do, largely because of this difference in the values 
of the coverin g fractions /*. ) If th e covering fraction is correlated with the accretion rate, 
as claimed by Batalha et al. ( 2002h . then the considerably larger accretion rates we measure 
compared to IMuzerolle et al.l (120001 ) should be associated with a larger covering fraction. It 
is interesting to note that the values of the energy flux and covering fraction we derive for 
TW Hya are very similar to those found by ICalvet fc Gullbringi (119981 ) for CY Tau, the star 
whose He I 1.083 fim line profile resembles that in our spectrum of TW Hya most closely. 

Despite the good agreement between our te mperature and density estimates and the 
model pred i ction s of ICalvet fc Gullbringi ( 119981 ) . our values do not agree with those of 
Bary et al.l (120081 ) for their sample of CTTSs, and our analysis rules out their values for 
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TW Hya with high confidence. iBary et all ( 120081 ) do not present results for individual stars, 
so it is unclear how much variation they find from object to object. The mass accretion rates 
listed for the objects in their sample are substantially higher than that measured for TW 
Hya. (The minimum value listed in their Table is more than double the value derived for 
TW Hya in §3.2.1 from the Pa (5 line and five times larger than that derived from the Br 7 
line.) However, one would expect that higher mass loss rates would generate larger volume 
densities. At the moment, we have no satisfactory explanation for this discrepancy. 

The high densities we derive provide a reasonable explanation for why other lines typ- 
ical of ionized regions (e.g., forbidden lines of O II and O III) are weak or absent in op- 
tical and NIR spectra of TW Hya obtained with narrow slits centered on the star (e.g., 
Rucinski &: Krautterl Il983l ). Since these ionization states are found primarily within the 
region where H is fully ionized and the density is well above the critical density for these 
transitions, the ground-state levels giving rise to these transitions are collisionally de-excited. 
Spectra obtained with wide slits, encompassing lower density r egions farther from the star, 
might be expected to reveal the usual forbidden lines (see e.g., iNajita et al.ll2010l ). 



Our analysis suggests that the O I lines seen in the NIR spectrum are produced by Ly 
(3 fluorescence. It is not immedia t ely ob vious, however, where the O I emission is generated. 
Muzerolle. Hartmann. &: Calvetl ( Il998af ) demonstrated that the profile of the 0.8446 line 
correlates with that of Pa 11 in a small sample of CTTSs, and therefore suggested that the 
lines were formed in a common region. Furthermore, the O I lines in our spectrum of TW 
Hya are unresolved at our resolution, which indicates that they are not arising from a region 
experiencing substantial infall or outflow, and therefore it is tempting to suggest that they 
are formed at the base of the accretion column, along with the H lines. However, production 
of the O I lines by Ly j3 fluorescence requires a high Ly /3 flux density in a region that contains 
neutral oxygen. The latter require ment is inconsistent with co mpletely ionized H gas in the 
pre-shock region in the models of ICalvet fc Gullbringi (119981 ) . A partially ionized region 
at the outer edge of the pre-shock region could be responsible for the O I emission. This 
region could also give rise to the [O I] 6300, 63 63 lines that are seen in the optical spectrum 
( lAlencar &: Batalhal 120021 ; iHerczeg et al.l 120071 ). However, fluorescenc e also requires a large 
optical depth in the Balmer transitions. Using the equati on given by iGrandil (Il980h for the 
escape probability of Ha, e^a and that by lCapriottil ( 119651 ) for the optical depth th q in terms 
of en Q , we derived a rough estimate of the Ha optical depth of thc ~ 84. Such a large optical 
depth in Ha seems incompatible with the finding that the observed H line st rength ratios are 
well reproduced by the theoretical Case B values (see also lBary et al.ll2008l ). which assumes 
that all transitions above the Ly series are optically thin. A more detailed model of the 
emission line region and its spectrum will be needed to resolve this puzzle. 
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4.4. Excess Spectrum 

Finally, our derived NIR excess flux spectrum, usually interpreted as arising from the hot 
inner accretion disk, does not match that from any recent model of TW Hya. While the H- 
band peak in our excess flux spectrum suggests a dust temperature of ~ 1600K, which would 
be expected for dust in the inner disk, the sharp drops in flux at both shorter and longer 
wavelengths are not reproduced by current optically thin disk models. In fact, as discussed 
above, we suspect that most of this excess emission does not arise from the inner accretion 
disk, but is rather due to the lower H~ opacity in the star's atmosphere relative to that of an 
older, higher gravity field dwarf. Therefore, during our observations, any NIR emission from 
the hot inner accretion disk must have been very small, much smaller than either previous 
estimates or the predictions of current models. It is possible that previous values of the NIR 
excess have been overestimated as a result of both the adoption of an incorrect spectral type 
for TW Hya, as well as the spectral differences due to differences in the opacity, gravity, and 
age between TW Hya and any comparison field dwarf. Matching a library template with a 
spectral type earlier than M2V to the short wavelength end of the NIR spectrum would yield 
substantially more excess emission at wavelengths > 2 /im than we estimate. In addition, 
scaling earlier spectral type templates to match the observed TW Hya flux levels between 
0.95 and 1.34 /xm yields an excess spectru m for which synthe tic photometry is in reasonable 



agreement with some of the models (e.g. Eisner et al.l 120061 ) . Nevertheless, in these cases 



the resulting excess flux spectra appear to be unphysical (because of negative values), and 
are certainly not the smooth blackbody functions characterized b y a single dust temper ature 



expected from the models or observed in other CTTS sources (e.g. lMuzerolle et al.ll2003l ). We 
caution that any attempts to determine the NIR excess emission arising from accretion disks 
around young stellar objects similar to TW Hya must account for any possible differences in 
age and gravity between the target objects and the spectral templates. Failure to do so could 
lead to spurious results. As we have shown here, the use of an old, higher gravity field dwarf 
as a spectral template can lead to an excess spectrum that primarily reflects differences in 
the H~ atmospheric opacity between the CTTS and the template, not the emission from an 
accretion disk. 



5. Conclusions 

We have presented new high signal-to-noise, moderate resolution NIR spectra of TW 
Hya obtained with SpeX at the IRTF. These data have enabled us to revise the spectral 
classification of this object. We find that the spectrum is more consistent with a spectral 
type of M2.5V, rather than the commonly adopted K7V. Based on this spectral classification, 
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we derive new values for the effective temperature (3400 K), mass (O.4M ), radius (1.29_R Q ), 
and age (3 Myr) of the system. Subtraction of an M2.5V template drawn from the IRTF 
Cool Star Spectral Library reveals emission lines of H, He, Ca, and O. The H lines consist of a 
narrow component, a much weaker broad component, and a red absorption feature. Analysis 
of the strengths of the narrow H lines indicates that they arise from a hot (T e ~ 20000 K), 
extremely dense (n e ~ 10 13 cm -3 ) region with a filling factor of /* > 0.02 and a thickness of 
AR ~ 10 2 — 10 3 km. The temperature and thickness values agree rema rkably well with the 



predi ctions from the model of magnetospheric accretion developed by ICalvet fc Gullbring 



( 119981 ). We interpret the structure in the line profile as a reflection of the relative locations 
of the emitting regions: the narrow core arises from the base of the accretion shock column 
at the stellar photosphere, while the inverse P Cygni absorption arises from the infall of gas 
in the accretion column and the extended wings arise from the outflow of material in a disk 
or coronal wind. Analysis of the O I permitted lines indicates that they are generated by Ly 
fluorescence. The location of the mostly neutral region generating the O I lines, however, 
is uncertain. We also find very little excess continuum emission arising from the optically 
thin hot inner accretion disk. The spectrum of the exc ess emission we der ive does not appear 



to agree with the model of the inner disk presented by lEisner et al.l ( 120061 ). Rather the shape 
of the excess suggests that it results from the lower H~ opacity in the atmosphere of TW 
Hya, due to the lower logg and younger age, compared to the older, higher logg value for 
field M2.5V stars used for the spectral comparison. 

WDV thanks M. C. Cushing for providing software useful for some of the analysis and 
plots presented in this paper, as well as for discussions about M stars. We also thank J. 
Najita for discussions regarding TW Hya, C. Johns-Krull for providing H a equivalent width 
measurements, and C. Howard for comments on the manuscript. 
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Wavelength ((im) 



Fig. 1. — Flux-calibrated NIR spectrum of TW Hya obtained with SpeX on the IRTF. 
The strongest emission lines are identified. Regions of poor atmospheric transmission 
(< 20%) are shown in dark gray, while regions of moderate atmospheric transmission 
( < 80%) are shown in light gray. This figure can be directly compared with the figures 



in 



Rayner. Cushing fc Vaccal (120091 ). The bottom plot shows the statistical S/N. 
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Fig. 2a. — (Top) Comparison of the normalized 0.8 — 1.1 fim spectrum of TW Hya with 
those of late K and early M dwarfs stars in the IRTF Library. Note the variation in the 
strength of the FeH absorption with spectral type. The identification of o ther f eatures can 
be found in iRayner. Cushing &: Vaccal (120091 ); ICushing. Rayner. fc Vaccal (120051 ). (Bottom) 
Residual spectrum obtained by subtracting the scaled M2.5V spectrum from that of TW 
Hya. 
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Fig. 2b. — (Top) Comparison of the normalized 1.1 — 1.35 /im spectrum of TW Hya with 
those of late K and early M dwarfs stars in the IRTF Library. Note the variation in the 
strength of the K I and H 2 absorption features with spectra l type. The identification of 
other features can be found in lRayner. Cushing fc Vaccal (120091 ); ICushing. Rayner. &: Vacca 
(120051 ). (Bottom) Residual spectrum obtained by subtracting the scaled M2.5V spectrum 
from that of TW Hya. 
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Fig. 2c. — Comparison of the normalized 1.4 — 1.8 /im spectrum of TW Hya with those 
of late K and early M dwarfs stars in the IRTF Library. Note the variation in the strength 
of the Mg and Al absorption features, as well as the complex of Fe and Si lines between 
1.60 and 1.65 Aim, with spectral type. The identification of other features can be found in 
Rayner. Cushing &: Vaccal ( 2009 ). (Bottom) Residual spectrum obtained by subtracting the 
scaled M2.5V spectrum from that of TW Hya. 
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Fig. 2d. — Comparison of the normalized 1.95 — 2.5 /im spectrum of TW Hya with those 
of late K and early M dwarfs stars in the IRTF Library. Note the variation in the strength 
of the Mg absorption feature with sp e ctral t ype. The identification of other f eatures can 
be found in iRayner. Cushing &: Vaccal (120091 ) ; ICushing. Rayner. & Vacca J2005h . (Bottom) 
Residual spectrum obtained by subtracting the scaled M2.5V spectrum from that of TW 
Hya. 



-33 - 




K I 1.253 |jm 



K5 MO M5 LO 



15 



10 



•< 

LU 



Na I 1.140 urn 



1 



5g 



5 



KO K5 MO M5 LO 



6 i 5 



3 4 

LU 



2 - 



Mg I 1 .504 pirn 



Mi 



KO K5 MO M5 LO 




KO K5 MO M5 

Spectral Type 



>' 1 V 



-t- - - - -A 



5 2? 
UJ : J 



5} 



1 : 



Mg I 1.711 urn 



KO K5 MO M5 

Spectral Type 



LO 



Fig. 3. — The strength of EW of various diagnostic spectral features as a function of spectral 
type for K-M dwarfs in the IRTF Spectral Library. The dashed lines denote the ±lcr range 
in the EWs for TW Hya. 
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Fig. 4.— Ratio of the EW of the Na I 1.140 fim feature to the EW of the Mg I 1.183 /im 
feature for K-M dwarfs from the IRTF spectral library. The dashed lines denote the ± la- 
range in the EW ratio for TW Hya. 
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Fig. 5. — Contours of x 2 values (relative to Xmin) as a function of electron density n e 
and temperature T e from fits of the theoretical H line fluxes to the values derived from the 
continuum subtracted observations of TW Hya. Contour levels are 1.05, 1.1, 1.2, 1.3, 1.5, 2.0, 
and 3.0 Xmin- Pluses denote the temperature and density combinations for which theoretical 
line ratios were computed and compared to the observations. The solid triangle denotes the 
best fit. 



-36 - 



1.00 



x 0.10 



CD 

sz 



CD 
> 

CD 
DC 



0.01 



Pa * 








Rr ^\ ~ ~ ~ ~ - 
i 1 1 \ N ^ 

\^ ^^^^^ 

V\ ^\ 








1 P f <k 

$ 


J 

$ 

□ 




T e = 20000 K, N e = 10 13 crrf 3 

1 -- T e = 3000 K, N e = 1 10 cm" 3 


Hu , 

< 


> 

i 

> 





10 



12 



14 



N 



upper 



Fig. 6. — Observed (points) and best-fitting theoretical line fluxes (solid lines, corresponding 
to T e = 2 x 10 4 K and n e = 10 13 cm 3 ) relative to Pa /3. We detected only a single line of the 
Humphreys series (Hu 11); the solid green square denotes the theoretical value for this line. 
The dotted lines and o pen green s q uare d enote the theoretical line fluxes for the temperature 
and density found by iBary et al.l (120081 ) in their study of 15 CTTS systems (T e = 3000 K 
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Fig. 7. — Pa (3 (left) and Br 7 (right) line profiles. Note the strong core centered at zero 
velocity, atop a weaker and broader component, and a weak absorption component in the 
red wing at 300 km s -1 . The dotted line is a Gaussian with the FWHM given by the nominal 
spectral resolution. 
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Fig. 8. — He I 1.083 /zm line profile. Note the red and blue absorption components. The red 
absorption feature is located at the same velocity as the similar feature in the H profiles. 
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Fig. 9. — (Top) Comparison between the spectrum of TW Hya and that of HD 237903 (K7V; 
left) and Gl 381 (M2.5V, right). The K7V and M2.5V spectra have been scaled to match 
the flux levels of TW Hya between 0.9 and 1.35 /mi. (Middle) Excess flux spectra generated 
by subtracting the scaled K7V spectrum (left) and the scaled M2.5V spectrum (right) from 
the observed sp ectrum of TW Hya . The red dashed line denotes the excess flux predicted 
by the model of Eisner et al.l (120061 1 . (Bottom) The percentage excess relative to the scaled 
template spectra. 
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Fig. 10. — The excess flux spectrum (black) generated by subtracting a scaled M2.5V spec- 
trum from the spectrum of TW Hya, compared with the excess flux spectrum (red) generated 
by subtracting a scaled M2.5V spectrum from an M2.5III spectrum. The latter excess spec- 
trum has been scaled to match the TW Hya excess spectrum between 1.3 and 2.2 /im. 
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Fig. 11. — The observed NIR spectrum of TW Hya (black line, offset by 0.2) compared to 
a composite spectrum (grey) generated from the sum of an M2.5V (red) and a K7V (green) 
spectrum. 
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Fig. 12. — H-R Diagram for TW Hya (thick cross) . Overplotted are t he theoretical 
tracks (solid lines) an d isoc hrones (dash ed lines) of Baraffe et al.l ( 1998 ) (upper left), 
D'Antona fe Mazzitellil (119981) flower left). ISiess. Dufour. &: Forestini (upper right), 



and iTognelli. Prada Moroni. Sz Degl'Innocenti feoiol ) (lower right). 
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Table 1. EWs of Selected Absorption Features in the NIR spectrum of TW Hya 



Feature 


A 


EW 


a 




jLtm 


(A) 


(A) 


All 


1.313 


2.38 


0.08 


All 


1.675 


4.13 


0.17 


All 


2.110 


0.80 


0.08 


All 


2.117 


1.12 


0.08 


KI 


1.253 


0.72 


0.04 


Mg I 


1.183 


0.97 


0.03 


Mg I 


1.504 


4.81 


0.14 


Mg I 


1.488 


1.51 


0.13 


Mgl 


1.711 


2.61 


0.21 


Mgl 


2.107 


0.50 


0.07 


Nal 


1.140 


2.94 


0.34 


Nal 


2.206 


4.41 


0.08 


Nal 


2.336 


2.12 


0.16 


FeH 


0.99 


3.10 


0.08 
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Table 2. Fluxes of Emission Lines in the NIR spectrum of TW Hya 



A Q bs a Flux a Identification 
(jum) (10~ 17 W in" 2 ) 



0. 


.8183 


8.21 


0. 


,63 


Nal 


0. 


.8195 


11.18 


0. 


.71 


Nal 


0. 


,8448 


13.61 


0. 


.58 


OI 


0. 


.8502 


17.03 


0. 


.83 


Ca II + Pa 


0. 


,8545 


19.16 


0. 


.70 


Ca II + Pa 


0. 


.8602 


6.01 


0. 


.82 


Pa 14 


0. 


,8630 


1.85 


0. 


.56 


? 


0. 


.8665 


18.52 


0. 


.75 


Ca II + Pa 


0. 


.8753 


12.02 


0. 


.84 


Pa 12 


0. 


.8864 


16.83 


0. 


.80 


Pa 11 


0. 


.9018 


13.99 


0. 


.79 


Pa 10 


0. 


.9232 


32.16 


0. 


.88 


Pa 9 


0. 


.9549 


47.15 


1. 


.09 


Pa 8 (e) 


1. 


.0051 


64.00 


1. 


00 


Pa 7 (5) 


1 


.0406 


9.01 


1 


.04 


[N I]? 


1 


.0941 


92.11 


0. 


.92 


Pa 6 (7) 


1. 


.1289 


9.99 


0. 


.40 


I 


1. 


.2821 


176.89 


1. 


.79 


Pa 5 (/?) 


1. 


.6812 


2.45 


0. 


.87 


Br 11 


1. 


.7368 


8.69 


1 


33 


Br 10 


1. 


.9449 


10.11 


1 


.11 


Br 8 (e) 


2. 


.1661 


22.32 


1. 


.37 


Br 7 (7) 


3. 


.0395 


3.27 


0. 


32 


Pf 10 (e) 


3. 


.2971 


4.96 


0. 


.30 


Pf 9 (5) 


3. 


.7404 


7.70 


0. 


.47 


Pf 8 (7) 


4. 


.0523 


62.38 


0. 


.55 


Br 5 (a) 


4. 


.6538 


12.57 


0. 


.84 


Pf 7 03) 


4. 


.6727 


1.46 


0. 


.58 


Hu 11 (e) 



-44 - 



a Observed vacuum wavelength 



